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Abstract 

In this paper, we propose a new comparison tool for spatial homogeneity 
of point processes, based on the joint examination of void probabilities and 
factorial moment measures. We prove that determinantal and permanental 
processes, as well as, more generally, negatively and positively associated point 
processes are comparable in this sense to the Poisson point process of the same 
mean measure. We provide some motivating results and preview further ones, 
showing that the new tool is relevant in the study of macroscopic, percolative 
properties of point processes. This new comparison is also implied by the 
directionally convex (dcx) ordering of point processes, which has already been 
shown to be relevant to comparison of spatial homogeneity of point processes. 
For this latter ordering, using a notion of lattice perturbation, we provide a 
large monotone spectrum of comparable point processes, ranging from periodic 
grids to Cox processes, and encompassing Poisson point process as well. They 
are intended to serve as a platform for further theoretical and numerical studies 
of clustering, as well as simple models of random point patterns to be used 
in applications where neither complete regularity nor the total independence 
property are not realistic assumptions. 
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1. Introduction 

Usual statistical approach to the study of clustering in point processes (pp) consists 
in the evaluation of Ripley's K function, pair-correlation function, or contact distri- 
bution function (also called the empty space function). However, such a comparison 
of local characteristics seems a weak tool for the study of the impact of clustering 
on some macroscopic properties of pp such as those required in continuum percolation 
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models. We are particularly motivated by heuristics indicating that pp exhibiting more 
clustering should have larger critical radius for the percolation of its spherical-grain 
Boolean model than a spatially homogeneous pp. 

It was observed in [4], that the directionally convex (dcx) order on pp implies the 
ordering of K functions as well as pair-correlation functions, in the sense that pp larger 
in the dcx order have larger K functions and pair correlation functions, while having the 
same mean number of points in any given set. Unfortunately, the examples from [4] are 
mostly only some doubly-stochastic Poisson pp, which are dcx larger than Poisson pp 
(we call them super-Poisson in this article) . In order to provide more examples of dcx 
ordered pp, in particular smaller than Poisson (we call them sub-Poisson) , we study in 
this paper a notion of perturbation of a pp consisting of independent replication and 
translation of points from some given, original pp. A key observation is that such a 
perturbation is dcx monotone with respect to the convex order on the number of point 
replications. In particular, perturbing a deterministic lattice in the above sense, one 
can obtain examples of both sub- and super-Poisson pp, with the Poisson pp itself 
obtained when the number of point replications has a Poisson distribution. We believe 
these examples can be useful for modeling of real phenomena for which neither lattice 
nor Poisson assumptions can be justified. In this paper, we will also use them to 
illustrate the aforementioned heuristic on the impact of clustering on the percolation 
of Boolean models. 

However, many examples of pp considered as clustering less or more than the Poisson 
pp of the same intensity escape from the dcx comparison; For example, determinantal 
and permanental pp (cf. [3]). In fact, despite some structural similarities of these pp 
to the perturbed lattices, we are able to show for them dcx order only on mutually 
disjoint simultaneously observable sets, and not on all bounded Borel sets, required for 
the full dcx order. 

The properties of positive and negative association (cf [7,18]) are also used to define 
classes of pp that, respectively, cluster more or less than the completely independent 
(i.e., Poisson) pp. But it is not known if these properties imply or are implied by 
the dcx ordering with respect to Poisson pp. Though one suspects many pp such as 
determinantal or hard-core pp should be negatively associated, it is not known if they 
actually are 1 . 

In order to unify the approach to matter in hand and provide more examples of pp 
comparable to Poisson pp, we define two more classes of pp: weakly sub-Poisson — as 
pp having both void probabilities and factorial moment measures smaller than Poisson 



1 However, there are examples of negatively associated discrete measures including determinantal 
ones (see [12, Theorem 6.5]). 
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pp of the same mean measure, and weakly super-Poisson — as having these character- 
istics larger than Poisson pp of the same mean measure. It is almost straightforward to 
see that this new classification is indeed weaker than sub- and super-poissonianity based 
on the dcx ordering. We prove also that it is also weaker than association: positive 
association implies the weak super-poissonianity, while negative association implies the 
weak sub-poissonianity. A good news is that permanental and determinantal pp can 
be proved to be weakly super- and sub-Poisson respectively. Also, as it turns out many 
of the results can be proven under these weaker assumptions of weakly sub-Poisson or 
super-Poisson than association or dcx ordering. 

Paper organization The necessary notions, notations and basic facts are introduced 
and recalled in Section 2. In Section 3, we define classes of strongly and weakly 
sub- and super-Poisson pp and, as a main result, we prove that weak sub- or super- 
poissonianity is implied by negative or positive association, respectively. We study the 
perturbed-lattice pp in Section 4 and determinantal and permanental pp in Section 5. 
In Section 6, we discuss some further theoretical implications (especially percolation) of 
the presented ideas, and the modeling application. Lemma A.l, which is of independent 
interest and used in this paper for showing dcx ordering of perturbed lattices and 
determinantal and permanental point processes (on mutually disjoint simultaneously 
observable sets) is proved in the Appendix. 

2. Notions, notation and basic facts 

Point processes We assume the usual framework for random measures and point 
processes on d-dimensional Euclidean space M d (d> 1), where these are considered as 
random elements on the space M(IR ci ) of non-negative Radon measures on M. d (cf [11]). 
A point process (pp) $ is simple if a.s. <!>({:e}) < 1 for all x € M. d . We denote 
by v(B) = P($nB = 0), the void probabilities of pp $ and by a^(-), the factorial 
moment measure of Recall that for simple pp, a^ k \B\ x . . . x Bf.) — E Yii=i f° r 
pairwise disjoint bBs Bi (i = 1, . . . , k). The k th joint intensity, p^ : (R d ) k — > [0, oo) 
is the density (if it exists) of a^(-) with respect to the Lebesgue measure dxi . . . dx/.. 
Recall that the joint intensities p^ k \ k > 1 characterize the distribution of a pp. The 
above facts remain true even when the densities p( k > are considered with respect to 
IliLi /i(dxfc) for an arbitrary Radon measure p on R d . As always, a pp or a random 
measure on M. d is said to be stationary if its distribution is invariant with respect to 
translation by vectors in M. d . 

Directionally convex ordering A Lebesgue-measurable function / : R k —> R is 
said to be directionally convex (dcx) if for every x £ ]R fc , e, S > 0, i,j € {1, . . . , k}, we 
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have that A\A J s f(x) > 0, where A\f(x) :— f(x + eai) — f(x) is the discrete differential 
operator, with {ei}i<i<k denoting the canonical basis vectors for M. k . We abbreviate 
increasing and dcx by idcx and decreasing and dcx by ddcx (see [17, Chapter 3]). For 
real-valued random vectors of the same dimension X and Y, X is said to be less than 
Y in dcx order (denoted X < dcx Y) if E(f(X)) < E(/(F)) for all / dcx such that 
both the expectations are finite. For two pp on M. d , one says that $i(-) <dcx ( ■ ) ' if 
for any B fc bBs in R k , . . . , *i(B fc )) < dc:c (* 2 (Bi), ■ ■ ■ , Wfc)); cf [4]. 

The definition is similar for other orders, i.e., those defined by idcx, ddcx functions. It 
is enough to verify the above conditions for Bi mutually disjoint. In order to avoid 
technical difficulties, we will consider here only pp whose mean measures E($(-)) are 
Radon (finite on bounded sets). For such pp, dcx order is a transitive order 2 . 

It is easy to see that <dcx 3>2(0 implies the equality of their mean measures: 

E($i(-)) = E($2('))- Moreover, as shown in [4], higher-order moment measures are 
non-decreasing in dcx order on pp, provided they are cr-finite 3 . In addition, dcx 
ordering allows to compare also the void probabilities as stated in the following new 
result: 

Proposition 2.1. Denote by v\{-),V2{ - ) the void probabilities of pp $i and <E>2 on Mr 
respectively. If $1 <ddcx $2 then v\{B) < v 2 (B) for all bBs B C R d . 

Proof. This follows directly from the definition of dcx ordering of pp, expressing 
Vj{B) = E(/($ i (B))), j = 1,2, with the function f(x) = max(0, 1 - x) that is 
decreasing and convex (so ddcx in one dimension). □ 

In particular, the latter result implies ordering of all contact distribution functions 
(empty space functions) for pp comparable in dcx order and not having fixed atoms 4 . 
We see in the joint comparison of moment measures and void probabilities of pp having 
equal mean measures, a new tool for comparison of their clustering properties, weaker 
than dcx order but more easy to verify. 

Positive and negative association Denote by Cov (XY) — E(XY) — E(X) E(Y) 
covariance of random variables X, Y. A point process $ is called associated if 
Cov (/($(Bi), ... ,$(-B fe )), ... ,$(B k ))) > for any finite collection of bBs 

Bi, . . . , Bk C lR d and /, g continuous and increasing functions taking values in [0, 1]; 
cf [7] . This property is also called positive association, or the FKG property. The theory 



2 Due to the fact that each dcx function can be monotonically approximated by dcx functions fi(-) 
which satisfy fi(x) = ()( x . x ) at infinity, where H^Hoo is the Loo norm on the Euclidean space; 
cf. [17, Theorem 3.12.7]. 

3 <r-finiteness condition is missing in [4]; see [24, Prop. 4.2.4] for the correction 

Satisfying Pr{x £ *} = 1 for all x e R d . 



4 



for the opposite property is more tricky, cf [18], but one can call $ negatively associated 
if Cov . . . , $(£*.)), g($(B k+ i), $(Bi))) < for any finite collection of bBs 

Bi , . . . , Bi C R d such that (B x U • • • U B k ) H (B fc+1 U • • • U Bj) = and /, g increasing 
functions; Both definitions can be straightforwardly extended to random measures. 

3. Comparison of clustering to Poisson pp 

We call a pp sub-Poisson (respectively super- Poisson) if it is smaller (larger) in dcx 
order than the Poisson pp (necessarily of the same mean measure). (More precisely 
we should have called these processes dcx-sub-Poisson or cfcx-super-Poisson pp, but we 
omit the word dcx for simplicity.) Examples of such pp are given in Section 4. A weaker 
notion of sub- and super-poissonianity can be defined when comparing only moment 
measures or void probabilities. Bearing in mind that Poisson pp can be characterized 
as having void probabilities of the form v[B) = exp(— a(B)), where a(-) is its mean 
measure, we say that a pp <& is weakly sub-Poisson in the sense of void probabilities 
(j/-weakly sub-Poisson) if 

P = 0) < e -**t B » (1) 

for all Borel sets B C R d . Similarly, we say that a pp $ is weakly sub-Poisson in the 
sense of moment measures (a-weakly sub-Poisson) if 

k k 
aW(B i )<l[aU(B i ) = l[E(B i ) (2) 

4 = 1 1=1 

for all mutually disjoint bBs Bi C K d . When the inequalities in (1) and (2) are reversed, 
we will say that $ is i/-weakly super-Poisson or a-weakly super-Poisson respectively. 

Finally, we will say that $ is weakly sub-Poisson if $ is a-weakly sub-Poisson and v- 
weakly sub-Poisson. Similarly, we define weakly super-Poisson pp. Examples of weakly 
sub- and super-Poisson pp are given in Section 5. 

The fact that dcx ordering implies ordering of moment measures and void prob- 
abilities lend credence to our usage of the terms weak sub- and super-Poissonianity. 
Interestingly, these inequalities are also implied by negative and positive association. 
The following result is a key observation in this matter. 

Proposition 3.1. Consider point process $ with Radon mean measure a(-) = E($(-)). 
If <fr is simple, has Radon second-order factorial moment measure a^(-) and 

Pr{ - 0, $(B 2 ) = } < Pr{ $(Bi) = } Pr{ *(B a ) - }, (3) 

for any two disjoint bBs B\ and Bi, then $ is v -weakly sub-Poisson. 
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If the mean measure a(-) of $ is diffuse (without atoms) and <f> satisfies (3) with 
the reversed inequality (>) for any two disjoint bBs B\ and Bi, then $ is v -weakly 
super-Poisson. 

Proof. Define a set function Q(B) = — log(Pr{$(£>) = 0}). Regarding the first 
statement, it is immediate to see that Q is non-negative and, under assumption (3) 
super-additive; i.e., for any finite k > 1 and any pairwise disjoint bBs Bj, j = 1, . . . , k 
Q(Bi U . . . U Bk) > Xw=i Q (Bj)- In order to prove the result, we need to show that 
Q(B) > a(B), for any bBs B. To this regard note by the super-additivity of Q that 
for any bBs B 



where the "sup" is taken over all finite partitions of B into bBs Bi. Moreover, for any 



Pr{$(B) = 0} = 1 - E($(S)) + E($(B)1($(B) > 2)) - Pr{$(B) > 2} 
< 1-E($(B)) + E($(B)($(B)-1)+) 
= 1 - a(B) + a {2) (B x B) 

and hence Q(B) = - log(Pr{ = 0}) > a(B) ~ (B x B) . Consequently, by (4), 

for any bBs B 

Q(B) > supJ2 (a{Bj) - n (2) {B 3 x fl,-)) = a ( B ) - inf X) aW ^ x B ^ ' 

3 jeJ J jeJ 

due to finiteness of all terms. In order to complete the proof it is enough to show 
that the "inf" term is equal to zero. To this regard, for a given e > define = 
{B x B 3 (x,y) : \{x,y) — [z,z)\ < e for some z € B}. Note that can be seen as 
some neighborhood of the intersection of the diagonal with B x B. Note also that for 
any e > there exits a suitable fine partition / of B such that ^2j e j ct^ (Bj x Bj) < 
a( 2 \A e B ). (For example, take a finite coverage of B by balls of radius e, which exists 
by local-compactness of the space, and refine it to have disjoint partition of B.) By the 
local finiteness and a-additivity of a^ 2 \ lim e ^ (A e B ) — a^ 2 \{(z,z) : z g B}) = 0, 
where the last equality follows from the assumption that <& is simple. This completes 
the proof of the first statement. 

For the second statement, we will show that Q(B) < a(B). To this regard, note 
that the reversed inequality in (3) implies that Q(-) is sub-additive and consequently, 




ieJ 



(4) 



bBs B 
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for any bBs B, 

Q(B) =ini J2Q(Bj), (5) 

where "inf" is over all finite partitions of B. Moreover, observe that Pr{ Q(B) = } > 
l-a(B) and that, for < x < e, -log(l-x) < x(l + 5(e)), where (5(e) = e/(2(l-e) 2 , 
which can be shown by the Taylor expansion with Lagrange form of the remainder 
term of order 2. Since a(-) is diffuse, for any e > there exists a partition J of bBs B 
such that ct(Bj) < e for all j € J. For such a partition J, 

Q(B) <J2~ Ml " < a(B)(l + 5(e)) . 

je.J 

The proof follows from the observation that 8(e) — > when e 0. □ 

Corollary 3.1. ^4 negatively associated, simple pp with a Radon mean measure is 
weakly sub-Poisson. A (positively) associated pp with a Radon, diffuse mean measure 
is weakly super-Poisson. 

Proof. Inequality (2) or its inverse (i.e.; a- weak sub- or super-poissonianity) follows 
directly from negative association or association, respectively. The i/-weak sub- or 
super-poissonianity follows from Proposition 3.1. Indeed, inequality (3) or its inverse 
can be derived easily from negative association or association, respectively. Moreover, 
note by (2), that any factorial moment measure cr n '(-) of a simple, a-weakly sub- 
Poisson pp with Radon mean measure is also Radon. This completes the proof. □ 

In fact, negative association or (dcx) sub-Poissonianity implies something stronger 
than a-weak sub-Poissonianity. Namely, we have that, a^ k+l ^(-) < aW(-)a (0 (-) for 
integers k, I > 0. Positive association or (dcx) super-Poissonianity imply the reverse 
inequality. 

4. Perturbed lattices and point processes 

It was observed in [4] that Poisson-Poisson cluster pp, Levy based Cox pp, Ising-Poisson 
cluster pp are super-Poisson pp. In this section, we present more examples of pp, which 
are dcx comparable to Poisson pp. We begin with a general model of a perturbation 
of a pp and prove our key result on the dcx ordering of such pp. 

4.1. Perturbation operator 

Let $ be a pp on R d and Af (■,-), X(-,-) be two probability kernels from R d to 
non- negative integers Z + and M. d , respectively. Consider the following independently 
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marked version of the pp $, $ per * = {(X, Nx , Yx)}xg$ where given $: 

• N x , X £E 4> are independent, non- negative integer- valued random variables with 
distribution P (N x G • | $) = •), 

• Yj = (Fjx :i = l)2 r ..), J6$ are independent vectors of i.i.d. elements of 
R d , with Kjx's having the conditional distribution P (Y iX € • | $) = •), 

• the random elements AT^, Yj are independent for all X 6 <t>. 

Consider the following subset of M d 

|J + (6) 

xe$ i=i 

where the inner sum is interpreted as when Nx — 0. The set <j> pert can (and will) be 
considered as a pp on R d provided it is locally finite. In what follows, in accordance 
with our general assumption for this article, we will assume that the mean measure of 
<j,pert j g i oca Uy finite (Radon measure) 

/ n(x)X(x, B - x) a{dx) < oo, for all bBs B C R d , (7) 

where a(-) is the mean measure of the pp $ and n(x) = J^feLi kN~(x, {A;}) is the mean 
value of the distribution Af(x, •). 

The pp $ per * can be seen as independently replicating and translating points from 
the pp <£, with the number of replications of the point X G $ having distribution 
J\f(X, •) and the independent translations of these replicas from X by vectors having 
distribution X(X, •). For this reason, we call $P ert a perturbation of $ driven by the 
replication kernel Af and the translation kernel X . 

An important observation for us is that the operation of perturbation of $ is dcx 
monotone with respect to the replication kernel in the following sense. 

Proposition 4.1. Consider a pp $ with Radon mean measure a(-) and its two pertur- 
bations $^ ert j = 1,2 satisfying condition (7), having the same translation kernel X and 
possibly different replication kernels Nj, j — 1,2, respectively. If Mix,-) < cx N~(x,-) 
( convex ordering of the conditional distributions of the number of replicas ) for a-almost 
allxeR d , then $f rt < dcx $f rt . 

Proof. We will consider some particular coupling of the two perturbations $^ er * , j = 

1, 2. Given $ and Y x = (Y iX ■ i = 1, • • •) for each let $ iX = {j"=A x + Y ix}, 

where N 3 X has distribution Afj(X, •), j = 1, 2, respectively. Thus <j>^ er * = J2xe<s> ®jx> 
j = 1,2 are the two considered perturbations. Note that given $^ er * can be seen 
as independent superpositions of &jx for X e $. Hence, by [4, Proposition 3.2(4)] 
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(superposition preserves dcx order) and [17, Theorem 3.12.8] (weak and L\ convergence 
jointly preserve dcx order), it is enough to show that conditioned on $, $xx <dcx 
for every X £ $. In this regard, given consider X £ $ and let Bi, . . . , Bk be 
mutually disjoint bBs and / : M. k — > R, a cfca; function. Define a real valued function 
.g : Z -> R, as 



g(n) := E /(sgn(n) ^(l[r ix € B x - X], . . . , l[Y iJC eB t - X]) 



where sgn(n) = -A- for fi^O and sgn(O) = 0. By Lemma A.l, g(-) is a convex function 
on Z and by Lemma A. 2 it can be extended to a convex function g(-) on R. Moreover, 
E(§(i\&)|#) = e(s(JV£)|*) = e(/($ 3 -x(Bi),...,^x(B*))|*) for J = 1,2. Thus, 
the result follows from the assumption jVjp < cx N^- □ 

Remark 4.1. The above proof remains valid for an extension of the perturbation 
model in which the distribution X(X. •) of the translations Yix depends not only on 
the location of the point X £ $ but also on the entire configuration $; X(X, •) = 
X(X, <E>, •), provided condition (7) is replaced by finiteness of J M<J J Rd n(x)X(x, <j),B — 
x) C(d(x, 0)), where C(d(x, <j))) is the Campbell measure of 

4.2. Examples 

Perturbed Poisson pp Let $ be a (possibly inhomogeneous) Poisson pp of mean 
measure a(dx) on R d . Let N{x, ■) = z\ = 1(1 £ •) be the Dirac measure on Z+ 
concentrated at 1 for all x £ R d and assume an arbitrary translation kernel X satisfying 
a pert (A) — J Rd X(x,A — x) a{dx) < oo for all bBs A. Then by the displacement 
theorem for Poisson pp, <&P ert is also a Poisson pp with mean measure aP ert {dx). 
Assume any replication kernel A/^af, •), with mean number of replications n-z{x) — 
Y^k=\ kNs(x, {k}) = 1 for all x £ R d . Then, by the Jensen's inequality and Proposi- 
tion 4.1, one obtains a super-Poisson pp <J>^ er *. In the special case, when Af%(x, •) is 
the Poisson distribution with mean 1 for all x £ R d , &^ rt is a Poisson-Poisson cluster 
pp which is a special case of a Cox (doubly stochastic Poisson) pp with (random) 
intensity measure A (A) = X^e* X(x, A — x). The fact that it is super-Poisson was 
already observed in [4]. Note that for a general distribution of $, its perturbation <&^ rt 
is also a Cox pp of the intensity A given above. 



Perturbed lattice pp Assuming a deterministic lattice 4> (e.g. $ = Z d ) gives rise 
to the perturbed lattice pp of the type considered in [20] . Surprisingly enough, starting 
from such a 4>, one can also construct a Poisson pp and both super- and sub-Poisson 
perturbed pp. In this regard, assume for simplicity that $ = Z d , and the translation 
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kernel X(x, •) is uniform on the unit cube [0, l) d . Let Af(x, •) be the Poisson distribution 
with mean A (Poi(X)). It is easy to see that such a perturbation $P ert of the lattice 
% d gives rise to a homogeneous Poisson pp with intensity A. 

4.2.1. Sub-Poisson perturbed lattices. Assuming for A/i some distribution convexly [ex] 
smaller than Poi(X), one obtains a sub- Poisson perturbed lattice pp. Examples are 
hyper-geometric HGeo(n,m,k), m,k < n, km/n = A and binomial Bin(n, A/n), 
A < n distributions 5 , which can be ordered as follows: HGeo(n,m, Xn/m) < cx 
Bin(m, X/m) < cx Bin(r, X/r) < cx Poi(X) for A < m < min(r, n); cf. [22] 6 . Specifically, 
taking J\fi(x, •) to be Binomial Bin(n, X/n) for n > A, one obtains a cfcx monotone 
increasing family of sub-Poisson pp. Taking X = n = 1 (equivalent to Af(x, •) ~ £i), one 
obtains a simple perturbed lattice that is dcx smaller than the Poisson pp of intensity 1. 

4.2.2. Super- Poisson perturbed lattices. Assuming for Mi some distribution convexly 
larger than Poi(X), one obtains a super-Poisson perturbed lattice. Examples are 
negative binomial NBin(r,p) distribution with rp/(l — p) — X and geometric Geo(p) 
distribution with 1/p— 1 = A , which can be ordered in the following way: Poi(X) < cx 
NBin{r 2l X/{r 2 + X)) < cx NBin(r u A/(n + A)) < cx Geo(l/(l + A)) < cx ^ Aj Geo{p 3 ) 
with ri < r 2 , < Aj < 1, J2j = 1 anc i \IVi = A + 1, where the largest distribu- 
tion above is a mixture of geometric distributions having mean A; cf. [22]. Specifically 
taking M 2 (x, •) to be negative binomial NBin(n, X/ (n + A)) for n = 1, . . . one obtains 
a dcx monotone decreasing family of super-Poisson pp. Recall that NBin(r,p) is a 
mixture of Poi(x) with parameter x distributed as a gamma distribution with scale 
parameter pj (1 — p) and shape parameter r. 

From [15, Lemma 2.18], we know that any mixture of Poisson distributions having 
mean A is cx larger than Poi(X). Thus, the super-Poisson perturbed lattice with such 
a replication kernel (translation kernel being the uniform distribution) again gives rise 
to a Cox pp. 

4.2.3. Associated point processes: From [7, Th. 5.2], we know that any Poisson center 
cluster pp is (positively) associated. This is a generalization of our perturbation (6) 
of a Poisson pp $ (cf. Section 4.2) having form $ ci « s * e '' = J2xe<s>i X + with 
$x being arbitrary i.i.d. (cluster) point measures. Other examples of associated pp 



5 Bin(n,p) has probability mass function PBin(n,p) (*') = — p) n ~ l (* = 0, ...,n). 

HGeo(n, m, k) has probability mass function PHGeo(n,m,k) (*) = (™) Ct~™) / ( t) (max(fc — n + m, 0) < 
i < m). 

6 One shows the logarithmic concavity of the ratio of the respective probability mass functions, 
which implies increasing convex order and, consequently, cx provided the distributions have the same 
means. 

T PGeo(p)(i) = PC 1 -P)', PNBin(r,p)(i) = C" 1 ']' 1 )?' i 1 ~ VY ■ 
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given in [7] are Cox pp with intensity measures being associated. (It is easy to see by 
Jensen's inequality that all Cox pp are ^-weakly super-Poisson.) 

It is easy to see that the pp formed by throwing n i.i.d. points in a bounded 
region forms a negatively associated pp. Further, one can show that independent 
superposition of negatively associated pp is a negatively associated pp. Hence, simple 
perturbed lattices (cf. Section 4.2.1) are negatively associated. 

5. Determinantal and permanental point processes 

In the remaining part of this section, we focus on spatial determinantal and per- 
manental pp. We will show that they are, respectively, weakly sub- and super-Poisson 
pp. Some partial dcx comparison of these pp with respect to Poisson pp, namely on 
mutually disjoint, simultaneously observable sets, will be proved as well. 

5.1. Definition 

To make the paper more self-contained, we will recall a general framework from [3, 
Chapter 4], which allows us to study ordering of determinantal and permanental pp 
more explicitly; see also [2] for a quick introduction to these pp. 

Let K : M. d x R d — > C (where C are complex numbers) be a locally square-integrable 
kernel, with respect to /i® 2 on R 2d 8 . Then K defines an associated integral opera- 
tor JCd on L 2 (Z),/i) as ICafix) = j D K(x,y)f(y) fi(dy) for complex- valued, square- 
integrable / on D (/ € L 2 (D,fi)). This operator is compact and hence its spectrum 
is discrete. The only possible accumulation point is and every non-zero eigenvalue 
has finite multiplicity. Assume moreover that for each compact D the operator JCd is 
Hermitian 9 , positive semi-definite 10 , and trace-class; i.e., J^j \ < 00 j where Aj 5 
denote the eigenvalues of JCd . By the positive semi-definiteness of JCd , these eigenvalues 
are non-negative. 

Determinantal pp A simple pp on M. d is said to be a determinantal pp with a kernel 
K(x, y) with respect to a Radon measure [i on M d if the joint intensities of the pp with 
respect to the product measure \i® k satisfy p^(xi, . . . ,Xk) = det (K(x i: Xj)) 1< . . <k 
for all k, where (o'ij) 1<i - <k stands for a matrix with entries ay and det (•) denotes the 
determinant of the matrix. Note that the mean measure of the determinantal pp (if it 
exists) is equal to a(-) = J K(x,x) fi(dx). Assuming that the kernel K is an integral 
kernel satisfying the assumptions given in Section 5.1, the above equation defines the 



8 i.e., J D f D \K(x, y)\ 2 fi(dx) ^i(dy) < oo for every compact D C K d 

9 i.e., f D f(x) K. D g(x) fi(dx) = f D g(x)K. D f(x) /j,(dx) for all f,g 6 L 2 (D,fi) 

°i.e., f D f(x)IC D f(x)ii(dx) > 
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joint intensities. Then, there exists a unique pp $ de * on K d , such that for each compact 
D, the restriction of $ det to D is a determinantal pp with kernel Ko if and only if the 
eigenvalues of ICd are in [0,1]. 

Permanental pp Similar to the determinantal pp, one says that a simple pp is a per- 
manental pp with a kernel K(x, y) with respect to a Radon measure /i on M. d if the joint 
intensities of the pp with respect to [i® k satisfy p( k >(xi, . . . , Xk) — per(K(xi, Xj)) t< . . <k 
for all k, where per( • ) stands for the permanent of a matrix. Note that the mean 
measure of the permanental pp is also equal to a(-) = J K(x,x) fi(dx). Again, will 
assume that K(x,y) is an integral kernel. Then, there exists a unique pp $P erm on 
R d , such that for each compact D, the restriction of $P er ™ 1 to D is a permanental pp 
with kernel Kjj; cf. [3, Corollary 4.9.9]. We will call this pp permanental pp with the 
trace-class integral kernel K(x,y). From [2, Proposition 35 and Remark 36], we also 
know that <J>P erm is a Cox pp. 

5.2. Comparison results 

The following properties hold true for determinantal and permanental pp with a 
trace-class integral kernel K(x,y). 

Proposition 5.1. <J> det is a-weakly sub-Poisson, while $>P erm is a-weakly super-P oisson; 
both comparable with respect to the Poisson pp with mean measure a(-) given by a(D) = 
J D Kd(x,x) n(d) — J^j ^f> where the summation is taken over all the eigenvalues Aj 3 
ofK D . 

Proof. Since Kjj(x,y) is Hermitian and positive semi-definite, by Hadamard's 
inequality, det (K D (x i ,x j )) 1< . j<k < IjjLi K D (x u x t ) which implies (2). For $P erm , 
the proof follows from the permanent analogue of the Hadamard's inequality (see [14]). 
□ 

Proposition 5.2. $> det is v -weakly sub-Poisson, while $P erm %$ v-weakly super- Poisson. 

Proof. It is known that for each compact D, $ det (D) = £\ Bin(l, Xf) and $P e ""(D) = 
J2j Geo{ 1/(1 + Aj 3 )) where the summation is taken over all eigenvalues \® of K-d and 
Bin(l, Aj 3 ) are independent Bernoulli random variables while Geo(l/(l + Aj 3 )) are in- 
dependent geometric random variables; cf. [3, Theorems 4.5.3 and 4.9.4]. Consequently 

® det {D) < cx Poi(J2*f ) <cx ^ erm {D) , 

3 

with the left inequality holding provided $ det exists (i.e.; Aj 3 € [0, 1] for all compact 
D). Noting that convex order of integer- valued random variables implies ordering of 
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probabilities of taking value concludes the proof (see Proof of 2.1). □ 

Alternatively, one can prove the above result via Proposition 3.1 as the inequality 
(3) has been proved for determinantal pp in [9, Cor. 3.3.]. 

Corollary 5.1. Combining results of Proposition 5.1 and 5.2 we conclude that $ det is 
weakly sub-Poisson, while QP erm is weakly super-Poisson. 

In the next result, we will strengthen the above corollary, proving dcx ordering of 
finite-dimensional distributions of $ de * and <j)P erm on mutually disjoint simultaneously 
observable sets D±, . . . ,D k . Simultaneous observability means that the eigenfunctions 
of /Cy£). , restricted to Di are also eigenfunctions of K,r> i for every i — 1, . . . ,k. 

Proposition 5.3. Let $ det and $ per be, respectively, the determinantal and perma- 
nental pp with a trace- class integral kernel K and with $ de * being defined only if the 
spectrum of JC^d is in [0,1]. Denote by <fr Pm the Poisson pp of mean measure a(-) 
given by ct(D) = y\ Aj 3 for all compact D, where the summation is taken over all 
eigenvalues Aj 3 of JCd ■ Let D\ , . . . , D k be mutually disjoint, simultaneously observable 
(with respect to the kernel K) compact subsets ofW 1 and D = [jDi. Then 

(^ det (D 1 ),...,^ et (D k )) < dcx (^ Poi (D 1 ),...,^ Poi (D k )) 

< dcx ^> per (D 1 ),...,^(D k )) . 

Proof. Let denote the eigenvalues of JCd { and Aj 3 = Y^i=i^f,ii j = 

1, . . . , J are the eigenvalues of JCd with J denoting the number of eigenvalues of 
JCd {J — oo and allowed and in the latter case the sum is understood as 0). 
From [3, Prop. 4.5.9], we know that (V et (L>i), . . . , $ de '(£> fc )) = V 7 : V; V , £ 
where Nj ~ Bin(l, Aj 3 ) and given iVj's, ^ i > 1 are independent multinomial vectors 
Mul{l,\f tl /\f,...,\fJ\f) n : It is easy to sec that (V oi (.Di), . . . , $ Poi (L> fe )) = 

E^iEt'iCj. where M j ~ Poi(Xf) with ^/s and XP/s as defined above. Due 
to the independence of £i,/'s and the assumption Ylj=i < 00 (local trace-class 
property of JCd), it is enough to prove for each j, that Yli=i <dcx YliJi This 
follows from the fact that Bin(l,p) < cx Poi{p) (see Section 4.2) and Lemmas A.l, 
A. 2. This completes the proof of the inequality for the determinantal pp. 

Regarding the permanental pp, 
Kj ~ Geo(l/(l + Aj 3 ) and are as defined above; see [3, Theorem 4.9.7]. Similar 

to the above proof, the required inequality follows from the ordering, 2j=i £'i j ^dcx 



11 Mul(k,pi, . . . , pfc) with < Pi < 1, i p,; = 1, has probability mass function 

PMul(n, P1 ,...,p k )( n i, •••>«*) = mi^'nfcl Pi 1 ■ • -Pfc fc for ni + . . . + n fc = n and otherwise. 
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YliJi^ij f° r au i; which follows from the fact that Poi(X) < cx Geo(l/(l + A) (see 
Section 4.2) and Lemmas A.l, A. 2. This completes the proof. □ 

Remark 5.1. The key observation used in the above proof was that the number 
of points in disjoint, simultaneously observable sets can be represented as a sum of 
independent vectors, which themselves are binomial (for determinantal) or Poisson 
(for Poisson) or geometric (for permanental) sums of some further independent vectors. 
This is exactly the same representation as for the perturbed pp of Section 4.1 (available 
for any disjoint sets); cf the proof of Proposition 4.1. In both cases, this representation 
and Lemmas A.l, A. 2 allow us to conclude dcx ordering of the corresponding vectors. 

Example of the Ginibre process. Let $ be the determinantal pp on M 2 with kernel 
K((x!,x 2 ), (2/1,2/2)) = exp[(xiyi + x 2 y 2 ) + i{x 2 y\ - X1IJ2)], ZjtUj € M, j = 1,2, with 
respect to the measure [i(d(xi,x 2 )) = ir^ 1 cxp[— x\ — x%\ dx\dx 2 . This process is known 
as the infinite Ginibre pp. Denote by ^ G — {\Xi\ 2 : Xi 6 $ }, the pp on K + of the 
squared radii of the points of <& G . We know that an arbitrary finite collection of the 
annuli centered at the origin D, = {{x\,x 2 ) : ri < x\ + x\ < R{\ is simultaneously 
observable for this pp; cf. [3, Example 4.5.8]. Using this observation, Proposition 5.3 
and the fact that dcx order of pp on R is generated by the semi-ring of intervals, we 
conclude that ^ G is dcx smaller than the Poisson pp <i> 1 of unit intensity on R + . A 
partial result, * G ([0,r]) < cx $!([(), r]) for all r > 0, was proved in [4]. 

6. Applications and further research 

In what follows will give some motivating results and preview further ones motivating 
the ideas presented in this paper. 

6.1. Continuum percolation 

The Boolean model on a pp $ with radius r is defined as C($, r) := Uxe$ Bx(r), 
where E>x{r) denotes the ball of radius r centred at X. By percolation, we mean 
the existence of an unbounded connected subset of the Boolean model. The critical 
radius for percolation is defined as r c ($) := inf{r : P(C($,r) percolates) > 0}. We 
mentioned in the Introduction a heuristic saying that clustering worsens percolation. 
Now, we can use some family of perturbed-lattice pp (cf. Section 4.2), monotone in dcx 
order, to illustrate this heuristic. Indeed, Figure 1 hints at ordering of the critical radii 
of dcx ordered pp. However, as shown in [5], this conjecture is not true in general: 
there exists a super- Poisson pp with the critical radius equal to 0. What was also 
shown there, is that weakly sub-Poisson pp exhibit a (uniformly) non-trivial phase 
transition in their continuum percolation model (i.e., admit uniformly non-degenerate 
lower and upper bounds for the critical radius) . Similar results regarding fc-percolation 
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Figure 1: Mean fractions 
of nodes in the two 
largest components of the 
Boolean models generated 
by perturbed-lattice pp 
with Bin(n,l/n) and 
NBin(n, 1/(1 + n)) as 
replication kernels, having 
fixed spherical grains of 
radius r. The replication 
kernels converge in n, from 
below and from above in 
dcx, respectively, to Poisson 
pp whose critical radius is 
depicted by the dashed line. 




and SINR-percolation models (arising in modeling of connectivity of wireless networks) 
hold for dcx-sub Poisson pp. 

In what follows, we will present some intuitions leading to the above results and 
motivating our special focus on moment measures and void probabilities in the previous 
sections. Specifically, we will introduce two newer critical radii r c ,f c , which act as 
lower and upper bounds for the usual critical radius: r c < r c <f c . We will show that 
clustering acts differently on these new radii: 

r c ($ 2 ) < r c ($i) < r c ($i) < r c (*i) < r c ($ 2 ) 

for $1 having smaller voids and moment measures than $2- This sandwich inequality 
tels us that <&i exhibits the usual phase transition < r c (Qx) < oo, provided $2 
satisfies a stronger condition < ^(^2) and r c ($2) < 00. Conjecturing that it holds 
for Poisson pp $2, one obtains the result on (uniformly) non-trivial phase transition 
for all weakly sub-Poisson $1 — the one proved in [5] in a slightly different way. 

6.1.1. Moment measures and percolation Let W m = [—m, m] d and define h m ^ '■ (IR. ) — > 
{0, 1} to be the indicator of the event that Xi, . . . , Xk G ($nW / m )' c , \xi\ < r, mf xG g\y m \x— 
Xk\ < f, \xi+i — Xi\ < r VI < i < (k — 1), where dW n denotes the boundary of set W n . 
Let 7V mj fc($, r) — Ylxi x k eR d hm,k{Xi, . . . , Xk) denote the number of distinct self- 
avoiding paths of length fc from the origin O € i 2 to the boundary of the box W m in 
the Boolean model and N m (Q, r) — J2k>i N m ,k(§, r) to be the total number of distinct 
self-avoiding paths to the boundary of the box. We define the following "lower" critical 
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radius: 

r c ($) := infjr : liminf E(iV m ($, r)) > 0} . 

m 

Note that r c ($) = inf{r : lim m P (JV m ($, r) > 1) > 0}, with the limit existing be- 
cause the events {7V m ($,r) > 1} form a decreasing sequence in m, and by Markov's 
inequality, we have that indeed r c (&) > r c ($) for a stationary pp $. 

Proposition 6.1. Let Cj — C(<$>j,r), j = 1,2 &e fwo Boolean models with simple pp 
of germs $>j, j = 1, 2, and a-finite k th moment measures <X- for all k > 1 respectively. 
If a[ k \-) < a^\-) for all k > I, then r c ($i) > r c ($ 2 )- -fri particular, for a stationary, 
a-weakly sub-Poisson pp $i of unit intensity we have that 9dr c ($>i) d > 1 where Qd is 
the volume of the unit ball. 12 

Proof. The proof relies on the following easy derivation of the closed form expres- 
sions for E(7V m ($„r)) ,j = 1,2: E(iV m (*,r)) = E fe >i E(JV m , fc (*,r)) and 

E(iVm jfe (<I>,r)) = / /i m ,fe(xi,...,a; fc )a^ ) ((ia;i,...,dxfc). 

J(R d ) fc 

For the second part of the proof note that the above summation over k can be taken 
over k > m r := [m/rj — 1, where \_a\ denotes the larges integer not larger than a. 
Indeed, the maximal distance that can be reached by a path of length k in C(Q,r) 
is (k + l)r and hence h mi k > 1 implies that k > m r . Consequently, for a-weakly 
sub-Poisson pp $ 

E(iV m ($,r)) < ^ / h mM {x 1 ,...,x k )dx l ...dx k 

k>m r 

where the second inequality follows by releasing the condition that x k is close to dW m . 
Thus, E(iV m ($,r)) < oo for 9 d r d < 1 and hence the result 9 d r c ($) d > 1. □ 

6.1.2. Void probabilities and percolation Though we are interested in the percolation 
of Boolean models (continuum percolation models), but as is the wont in the subject 
we shall use discrete percolation models as approximations. For r>0,i£ Mr, define 
the following subsets of K d : Q r :— (—j-, ^\ d and Q r (x) := x + Q r . We will consider 
the following discrete graph parametrized by n € N : L* d = (Z^,E* d ) is the usual 



12 Similar to open paths from the origin to dQ m , one can define an open path on the germs of <E> 
crossing the rectangle [0,m] X [0,3m] X ... X [0,3m] across the shortest side and define yet another 
critical radius r a (&) as the smallest r for which such a path exists with positive probability for an 
arbitrarily large m ( [16, (3.20)]). An analogous inequality holds true for this critical radius too. 
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close-packed lattice graph scaled down by the factor 1/n. It has = where Z is 
the set of integers, as the set of vertices and the set of edges E* d := {(zi,Zj) £ (ZfJ 2 : 
Qt(^)nQf(^)^0}. 

A contour in L* rf is a minimal collection of vertices such that any infinite path 
in L* d from the origin has to contain one of these vertices (the minimality condition 
implies that the removal of any vertex from the collection will lead to existence of an 
infinite path from the origin without any intersection with the remaining vertices in the 
collection). Let F„ be the set of all contours around the origin in L* d . For any subset 
of points 7 C M. d , in particular for paths 7 £ T^(K), T n , we define Q 1 — {J zGl Q n {z). 

With these notations, we can define the "upper" critical radius r c ($). 

r c = r c ($) := inf jr > : for all n > 1, ^ P (C($, r) n Q 7 = 0) < 00} . (8) 

7er„ 

It might be seen as the critical radius corresponding to the phase transition when the 
discrete model L* d = (Z^,E* d ), approximating C(<I>,r) with an arbitrary precision, 
starts percolating through the Peierls argument. As a consequence, r c ($) > r c ($) 
(see [6, Lemma 4.1]). The following ordering result follows immediately from the 
definition. 

Corollary 6.1. Let Cj — C{<5>j,r), j = 1,2 be two Boolean models with simple pp of 
germs j — 1, 2. If $1 has smaller voids probabilities than $2 then r c ($i) < r c {<&2)- 

6.2. Multiple coverage 

For a point process $, define the fc-covered set Cfe($, r) :— {x £ K d : ELXi, . . . , X m £ 
$,m > k 3 x £ f\_i Bxi(r)}- Heuristically, clustering should reduce the 1-covered 
region but increase the fc-covered region for large k and we present a more formal 
statement of the same. Expected volume of the fc-covered region is one of the important 
quantities of interest in sensor networks and our result has obvious implications regard- 
ing the choice of k or the point process in the context of sensor networks. We introduce 
another stochastic order to state the result. We say that two random variables X, Y 
are ordered in uniformly convex variable order (UCVO)(A < uv Y) if their respect 
density funtions /, g satisfy the following conditions : supp(f) C supp(g), f{-)/g{-) is 
an unimodal function but their respective distribution functions are not ordered i.e, 
F(-) ^ G(-) or vice-versa (see [22]) and where supp(.) denotes the support of a function. 
Denote by || A\\ the Lebesgue's measure of bBs A C R d . 

Proposition 6.2. Let<S>i and $2 be two simple, stationary pp such that <$>i(Bo{r)) < uv 
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$ 2 (-Bo( r ) f or r > 0. Then there exists k$ > 1 such that for any bBs W C R d 

V fe: i< fe < fco , E(||Cfc($i,r)nW||)>E(||C fc ($2 S r)nW||) and 
y k>ko , E(\\C k {$ lt r)nW\\) < E(||C fc ($ as r)nW||). 

Proof. Firstly note that E(||Cfe($j, r) fl W\\) = J w P (&i(B x (r)) > k) dx = 
|| W||P($i (B (r)) > k) for i = 1,2. Now it suffices to show that P ($i(-B (r)) > fc) - 
P (&2(Bo{ r )) > fc) changes sign exactly once in fc for k > 1. This is implied by the 
UCVO order (see [22, Section 2 and Theorem 1]). 

It is known that log-concavity of f / g implies UCVO order as well as convex ordering. 
We have used the latter implication in our examples for sub-Poisson (see Section 4.2.1) 
and super-Poisson perturbed lattices (see Section 4.2.2). We can take for $i any of 
the sub-Poisson perturbed lattices presented in this article or determinantal pp and 
$2 as a Poisson pp. We can also take $i to be a Poisson pp and $2 to be any of the 
super-Poisson perturbed lattices presented in this article or permenantal pp. 

6.3. Further applications 

6.3.1. Minimal spanning forest (MSF) In a recent work [10], the authors show the 
connectivity of some approximations of the MSF for the weakly sub-Poisson point 
processes (cf. the conjecture by Aldous and Steele from 1992 that the MSF of Poisson 
pp is almost surely connected, proved by Alexander for dimension d — 2 in 1995). 

6.3.2. Fist passage percolation Existence of arbitrarily large voids in Poisson pp was 
shown in [1] to be a reason of infinite end-to-end packet-delivery delays in a time- 
space SINR model, studied in the framework of first passage percolation problem. 
Superposing the Poisson pp with an independent lattice of arbitrarily small intensity 
makes the delays finite. The latter result remains true when using a simple perturbed 
lattice, in which case the superposition is an example of a (dcx) sub-Poisson pp. 
Generalization to an arbitrary sub-Poisson pp is an open question. 

6.3.3. Applications in modeling In the context of wireless networks, pp are used to 
model locations of emitters/receivers. An ubiquitous assumption when modeling base 
stations in cellular networks is to consider deterministic lattices (usually hexagonal). 
On the other hand, mobile users are usually modeled by a Poisson pp. Both the 
assumptions are too simplistic. In reality, patterns of base stations are neither perfectly 
periodic, due to various locational constants nor completely independent because of 
various interactions: social, human interactions typically introduce more clustering, 
while the medium access protocols implemented in mobile wireless devices (as e.g. 
CSMA used in the popular WiFi technology) tend to separate active users. One clearly 
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sees the interest in perturbed-lattice models in this context. We believe also that our 
work may lay the groundwork in other domains, e.g. in social and economic sciences, 
where one studies the impact of clustering on the macroscopic properties of models (cf. 
e.g. [8]). 

6.3.4. Further research Another motivation to study sub-Poisson perturbed lattices 
comes from their relations to zeros of Gaussian analytic functions (GAF), cf [21], whose 
points exhibit repulsion at smaller distances and independence over large distances. 
However, the points seem more regularly distributed than in Poisson pp ( [19]). This 
asks the question whether zeros of GAF are comparable in some sense to Poisson pp. 
Gibbsian pp is another well-known class of point processes, which depending on the 
nature of the potential would be more or less clustering. Super and sub-poissonianity 
(even in the weak sense) have not been studied yet for Gibbsian pp. Devising statistical 
tests for sub-Poissonianity would be desirable. 



Appendix 

The following result, similar to [15, Lemma 2.17] is used in the proof of Proposi- 
tions 4.1 and 5.3. 

Lemma A.l. Let ^ = . . . , £*) £ M. k , (i = 1,2, . . .) be independent, identically dis- 
tributed vectors of (possibly dependent) non-negative random variables. Suppose f is a 
dcx function on R k . Then, the function g defined on 7L by g(n) = E^f(sgn(n) X)'='i 
for n ^ and g(0) = is convex on Z. 

Proof. We will prove that g(n) has non- negative second differences 

gin, - 1) + g(n + 1) - 2g(n) > for all n £ Z (9) 

and use the first part of Lemma A. 2. To prove (9), define G(n,m) :— J2iL n +i £i f° r 
< n < m and G(n, n) := (0, 0, . . . , 0) € R k for n > 0. We have for n > 1, 



2g(n) - 2E(y(G(0,n 

= E(f(G(0,n-l) + G(n-l,n)))+E(f[G(0,n) 

= E(/(G(0,n-l) + G(n,n+l)) +E(/(G(0,n) 

= E(/(G(0,n-l) + G(n,n+l) +/(G(0,n))) 

< E(/(G(0,n-l))+/(G(0,n) + G(n,n + l))) 

= g[n-l)+g(n + l), 
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where for the third equality we have used mutual independence of G(0, n — 1), G(n — 
1, n), G(n, n+1) and the fact that G(n—1, n) and G(n., n+1) have the same distribution, 
while the inequality follows from the dcx property of / and the assumption £ 4 > 0. 
This proves (9) for n > 1. Similar reasoning allows to show (9) for n < —1. Finally, 
note that for n = 

2 5 (0) = 2/((0,...,0)) 

= e(/(-G(0,1) + G(0,1)) +/((0,...,0))) 

< E(/(-G(0,l))+/(G(0,l)] 
= 5 (-l)+. 9 (l), 

We will prove the following two technical results regarding convex functions. We 
were not able to find their proofs in the literature. 

Lemma A. 2. Let g{n) be a real valued function defined for all integer n G Z and 
satisfying condition (9). Then for all n > 2 

n n 

<?(X> fc i) <E A ^) ( l0 ) 

i=l i=l 

for all ki G Z and < Aj < 1, X)"=i ^< = 1 swc/i that X)"=i ^ ^- Moreover, 
function g(-) can be extended to a real valued convex function defined on real numbers K. 

Proof. As mentioned in [13, Section V.16.B.10.a] it is easy to see that (9) is equiv- 
alent to (10) with n — 2. Assume now that (10) holds true for some n > 2 (and all 
< A, < 1, ki G Z, i = 1, . . . ,n satisfying Yl7=i ^ = 1> S"=i e We will prove 
that it holds true for n + 1 as well. In this regard, define for a given k G Z and distinct 
(otherwise we use directly the inductive assumption) ki, ■ ■ ■ , fc n +i G Z, the following 
functions: 

, _ x ,x x x _ k - k n+l - K{h - fc„+l) 



k n fcn+1 
n-1 

A n +i = A n +i(Ai, . . . , A n _i) := 1 — Aj — A„(Ai, . . . , A„_i) 



i=i 

n-1 



F(Ai,...,A n _i) := J~] Afg(fcj) + A n (A x , . . . , X n -i)g(k n ) 
i=i 

+A„+i(Ai, . . . , A n _i)<?(fc n+ i) . 
Note that for any Ai, . . . , A n _i we have J2i=i — 1 an d J2i=i ^i^i — k - Consider the 



20 



following subset of the n — 1-dimensional unit cube C := < (Ai, . . . , A„_i) G [0, 1] 



< A„ < 1,0 < A„ + i < l|. The proof of the inductive step will be completed 
if we show that F(-) > g(k) on C . In this regard note that C is closed and convex. 
Assume moreover that C is not empty; otherwise the condition (10) is trivially satisfied. 
Note also that F(-) is an affine, real valued function defined on H. n_1 . Hence, by 
the maximum principle, the affine (hence convex) function — F attains its maximum 
relative to C on some point (A°, . . . , A°_ x ) G dC of the boundary of C. Consequently, 
we have F(-) > F(A°, . . . , A°_ x ) on C and the proof of the inductive step will be 
completed if we show that i^(A°, . . . , A°_ x ) > g(k). In this regard, denote A° = 
A„(A] ) , . . . , Xn-i) and A° +1 = A„ + i(A", . . . , A°_ x ). Using the continuity of the functions 
A„(-) and A„ + i(-) is not difficult to verify that (A°, . . . , A°_ x ) G dC implies A° = 
for some j = 1, . . . , n + 1. Thus, by our inductive assumption, F(\1, . . . , A°_j) = 
Y^i=ii=jLj ^id{ki) > g(k), which completes the proof of (10) for all n > 2. 

For the second statement, we recall the arguments used in [23] to show that a 
function satisfying (10) for all n > 2 (called globally convex function there) has a 
convex extension on K. In this regard, consider the epigraph epi(g) := {(k, [/,) G 
Z x R : [i > g(k)} of g and its convex envelope epi co (g). It is easy to see that 
epi co ( 5 ) = {(x,/i) G M. 2 : (i > YJU X i9(ki) for some h € Z, < Ai < 1, ^=1 A * = 1 
and Y^i=i ^iki = x }- Define g(x) :— inf{/x : (x, fi) G epi co (g)} for all x G M. The 
convexity of epi co ((7) implies that g is convex on E and the global convexity (10) of g 
implies that g is an extension of g. This completes the proof. 
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